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ABSTRACT 

Based on the numerical solution of the time-dependent relativistic Euler equations 
onto a fixed Schwarzschild background space-time, we estimate the accretion rate of 
radial flow toward the horizon of a test perfect fluid obeying an ideal gas equation of 
state. We explore the accretion rate in terms of the initial density of the fluid for various 
values of the inflow velocity in order to investigate whether or not sufficiently arbitrary 
initial conditions allow a steady state accretion process depending on the values of the 
pressure. We extrapolate our results to the case where the fluid corresponds to dark 
matter and the black hole is a supermassive black hole seed. Then we estimate the 
equation of state parameters that provide a steady state accretion process. We found 
that when the pressure of the dark matter is zero, the black hole's mass grows up to 
values that are orders of magnitude above 10 9 M Q during a lapse of lOGyr, whereas 
in the case of the accretion of the ideal gas dark matter with non zero pressure the 
accreted mass can be of the order of ~ IMo/lOGyr for black holes of 1O 6 M . This 
would imply that if dark matter near a supermassive black hole acquires an equation 
of state with non trivial pressure, the contribution of accreted dark matter to the 
supermassive black hole growth could be small, even though only radial accretion is 
considered. 
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1 INTRODUCTION 

The problem of formation and evolution of supermassive 
black holes (SMBHs) in the center of gigant elliptic and 
disk galaxies remains unsolved. The black hole growth 
is usually related to its coexistence with the surround- 
ing matter, both baryonic and dark matter. Some mod- 
els consider these black holes are the result of the evo- 
lution of seed black holes of various initial masses that 
grow through accretion, for instance, if there are seed black 
holes as the result of primordial gas clouds collapse at the 
early evolution of galaxies, the seeds would have masses 
of the order of 10 3 to 10 4 M o ( |e.g. Eisenstein fc Loeb 1995| 
|Koush iappas fc Bullock 2004[ ), whereas when the seeds are 
considered to be the result of the collapse of relativistic 
massive stars the seeds should have started with masses of 
the order of hundreds of solar masses ( |Heger et al. 2003] ); 
also combined models of the collapse of primordial gas and 
dark matter are assumed as black hole seeds of interme- 
diate masses (|Umeda et al. 2009j) . The study of bolomet- 
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ric quasar luminosity function at various distances indicates 
that the accreted mass of supermassive black holes should 
be baryonic (|Small fc Blandford 19921 |Hopkins et al. 2007| ), 
however no final answer about the amount of accreted 
dark matter by SMBHs is available, for instance in 
(Pcirani & Frcitas 2008) it is shown that dark matter con- 
tributes with at most 10% of the total accreted mass, which 
together with observations of the bolometric quasar luminos- 
ity confirms that baryons are the main matter component 
that feeds the black hole through accretion. 

Standard analyses of SMBHs growth are based on 
the study of the evolution of phase space distributions, 
for instance, if SMBHs are primarily fed by collisionless 
dark matter, or stars, the key problem to overcome is 
how to refill the loss cone ( |e.g. Lightman & Shapiro " 1977| 
I Zhao c t al. 2002). Previous work has shown that the 
timescale for this is too long for collisionless mat- 
ter to contribute significantly to black hole growth 
(Read & Gilmorc 20 03)) . though box orbits in a triax- 
ial potential might provide an interesting caveat to 
this argument (Hollwy-Bockclman & Sigurdsson 2006J) . An- 
other possibility is runway growth at the center of 
dense star clusters. This could produce intermediate 
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mass black holes, which would be seed black holes that 
grow through later gas accretion ([Ebisuzaki ct al. 20011 
Portcgies & McMillan 2002 I , which is expected to be fu- 
eled by mergers ( |e.g. Volonteri fc Rees 2005[ l . Concerning 
the non-radial accretion process, particles have to over- 
come the angular momentum barrier in order to get the gas 
to the very center of a galaxy (e.g. King & Pringie 2006). 
Mechanisms proposed include the "bars within bars" model 
UShlosman et al. 198 9 Bcgclman et al. 2006 ) and direct col- 
lapse of cold clouds ()Loeb fc Rasio 1994[) . 

On the other hand, non standard models consider the 
SMBHs grow through the direct accretion of self-interacting 
or collisional dark matter ([Ostriker 2000, ?). In this case, 
the self interaction cross section is tuned to ensure that the 
loss cone is refilled by diffusion on an interestingly short 
timescale. Collisional dark matter has also been studied at 
galactic scale as a possible solution to the galactic core 
problem ( |Spergel fc Steinhardt 2000[| and the co rrespond- 
ing gravitational collapse ( |e.g. Moore et al. 2000[ ). The col- 
lisional dark matter also finds its candidates from exten- 
sions to the standard model of particles, for example singlets 
de.g. Holz fc Zee 2001[ ). 

In this paper we present a totally different approach 
that considers both the collisional and the collisionless 
dark matter. It consists in the solution of the fully rel- 
ativistic time-dependent Euler equations for a fluid on a 
Schwarzschild black hole space-time, where the fluid plays 
the role of dark matter. We model the collisional and col- 
lisionless dark matter as a relativistic ideal gas with and 
without pressure respectively, in order to study the possibil- 
ity that the pressure itself suffices to stop the accretion and 
pull the system toward a stationary state that allows the 
SMBH to have the observed masses. We measure the accre- 
tion rate of dark matter into the black hole parametrized by 
the equation of state, initial dark matter density and initial 
infall velocity. 

Given that we only consider radial accretion, our re- 
sults for both the collisional and collisionless cases are up- 
per bounds to more general accretion situations where non- 
trivial impact parameters of the dark fluid particles are con- 
sidered, for instance those where rotation of dark matter 
particles is considered. 

The paper is organized as follows. In section [2] we de- 
fine the relativistic Euler equations describing the fluid on a 
Schwarzschild black hole, whereas in section[3]we explain the 
numerical methods we use; in section [4] we show our results 
for the extrapolation to the supermassive black hole plus 
dark matter system and in [5] we summarize our conclusions. 
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where p is the proper rest mass density, u M is the 4- velocity 
of the fluid and V M is the covariant derivative consistent 
with the 4- metric of the space-time l|Misner et al. 1973j> . 

The stress energy tensor representing a perfect fluid in 
General Relativity is given by 



phu^u" + pg u 



(3) 



where p, h and g^ v are respectively the pressure, the specific 
enthalpy and the 4-metric tensor of the space-time, respec- 
tively. 

We choose to write down the space-time metric in the 
3+1 decomposition fashion for a coordinate system (t,x l ), 
that is (jMisner et al. 1973)1 : 



ds 



-{a 2 - PiP^dt 2 + 2/3 i dx i dt + -yijdxUx^ (4) 



where the latin indices i, j label the three spatial coordi- 
nates of the space-time, a is called the lapse function and 
together with the shift vector j3 l determine the foliation of 
the space-time in space-like hypersurfaces; 7^ is the induced 
three metric that relates proper distances on the spatial hy- 
persurfaces |Misncr et al. 1973[1 . 

With this information, the unknown quantities of Euler 
equations are the following primitive variables: the rest mass 
density of the fluid p, its pressure p and its velocity measured 
by an eulerian observer v % defined in terms of the spatial part 
of the 4- velocity of the fluid by v 1 = fb- + 77 , where W is the 
Lorentz factor W = 1 . . Given that there are five 

unknowns and four equations the system is closed with an 
equation of state p = p(p, e), where e is the specific internal 
energy, defined by h — 1 + e+ p/p. 

Even though one can write down relativistic Euler equa- 
tions straightforwardly in terms of the primitive variables, 
we work with conservative variables that allow one to write 
down the relativistic Euler equations as a set of flux balance 
laws 



2 RELATIVISTIC EULER'S EQUATIONS 
2.1 In general 

In order to track the evolution of a perfect fluid on a non-flat 
space-time it is necessary to write down the general rela- 
tivistic Euler equations. For a generic fixed space-time these 
can be derived from the local conservation of the stress- 
energy tensor and the local conservation of the number 
of baryons: 



du dF(u) 

~di + dx i 
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where q is the vector of conservative variables, F' are the 
fluxes in the three spatial directions, S is a source vector. 
This approach is adequate to implement the finite volume 
methods described below. For the perfect fluid above, a 



flux balance system of equations is ( 
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where 7 = det('yy) is the determinant of the 3-metric and 
r%„ are the Christoffel symbols. 

2.2 Equations for the spherically symmetric case 

We now specialize to the case of a Schwarzschild black hole, 
for which the line element in Eddington-Finkelstein coordi- 
nates x M = [t, r, 9, (/>], is written for the line element (Q with 
the following gauge and 3-metric 



2M 
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where M is the mass of the black hole. 

Notice that the radial is the only non-trivial component 
of the shift vector j3 r and also the only non-zero velocity is 
v r . This implies that there are only three non-trivial balance 
law equations. In terms of these geometric quantities, and 
provided the space-time is spherically symmetric, the system 
of flux balance equations reads: 



D 

T 



which obey the balance law 
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Finally, we adopt an ideal gas equation of state, in which 
the pressure is given in terms of the rest mass density and 
specific internal energy by 



p = (r - l)pe. 



(12) 
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where F is the adiabatic index or the ratio of specific heats. 
3 NUMERICAL METHODS 

Solution of the equations. In order to solve the system of 
equations (|10[) and we provide initial data for the rest 
mass density and velocity of the fluid, and use the equation 
of state (|12p to calculate the initial pressure in case it is 
non zero. In all our production runs we use a constant den- 
sity profile as initial data and various values of the constant 
density and velocity. 

We evolve the initial data using a finite volume 
approximation of equations (j 101) together with High 
Resolution Shock Capturing methods involving approxi- 
mate Riemann solvers that are able to track the po- 
tential shocks formed during the evolution of relativis- 
tic fluid equations Qe.g. LeVeque 1992[ ). We specially use 
the HLL Riemann solver and a constant piece-wise re- 
construction of variables at the inter cell boundaries 
( |for a detailed d escripti on see Toro 1997| . 

We also implement an atmosphere that prevents the rest 
mass density to be small enough as to provoke that relativis- 
tic Euler equations diverge or develop unphysical results, be- 
cause when the density is extremely small the enthalpy di- 
verges which provokes the numerical implementation to fail, 
which combines with round off errors through arithmetic op- 
erations, cumulate during the evolution and finally produce 
approximations that are not accurate. The atmosphere we 
use is such that p — max(p, 10 -10 ), which is a value that 
allows the convergence of our numerical methods. The need 
of the atmosphere is in fact one of the reasons why it is not 
-at the moment- trivial to simulate evolutions of fluids with 
ultra low densities and the reason why our results on dark 
matter are actually extrapolations of numerical experiments 
with rather higher densities of matter. 

Boundaries and boundary conditions. We carry out our 
simulations on a finite domain < r exc ^ r r max . On the 
one hand, given that we use Eddington-Finkelstein coordi- 
nates, which penetrate the horizon, we apply the excision 
technique (jSeidel fc Sucn 1992J, that is, a chunk of the do- 
main well inside the black hole's horizon is removed from the 
numerical domain from < r ^ r exc in order to avoid the 
black hole's singularity and the steep gradients of metric 
functions near there; since the light cones inside the hori- 
zon point toward the singularity there is no need to impose 
boundary conditions at r = r exc , instead, the fluid simply 
gets out the domain toward the singularity through such 
boundary. We choose the excision radius to be at r exc = M , 
that is, half the Schwarzschild radius, which is both, far 
enough from the singularity at r = and provides a buffer 
zone M < r < 2M for the fluid inside the horizon to flow 
smoothly from the horizon toward the excised chunk of the 
domain. 

On the other hand, we implement an artificial outer 
boundary at a finite distance from the black hole at r — 
r m a X , and apply simple boundary conditions through the 
extrapolation of the conservative variables. We have found 
that such boundary condition makes the artificial boundary 
to behave appropriately, that is, we have verified that our 
calculations do not depend on the location of the external 
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boundary, which is the requirement that must be fulfilled by 
a boundary boundary condition. 

Diagnostics. In order to diagnose our simulations on 
the fly we implement detectors located at various radii in 
the numerical domain, that is, we define spheres where we 
calculate scalar quantities. In particular we track the accre- 
tion mass rate as a function of time. In order to do so we 
calculate the accretion mass rate specialized for spherical 
flow on a spherically symmetric space-time by 

Mace = -A<KT 2 pW (v r - Q , (13) 

at various spherical surfaces, including the black hole's event 
horizon. 

Since we are very interested in checking the validity of 
our numerical results we also monitor both, the conservative 
and the primitive variables of our system of equations in the 
whole numerical domain. Then we apply a self-convergence 
test of our variables in order to verify that our calculations 
carried out under a discretized domain converge to a cor- 
rect solution in the continuum limit. That is, we want to 
make sure that our implementation is consistent, which in 
practice means that when increasing the numerical resolu- 
tion, our numerical solution converges fast enough toward a 
continuous solution of the continuous equations, which is a 
highly non-trivial test ( |e.g. LeVeque 1 992 ) . 

Units. We use geometrical units G = c = 1, and 
time and spatial units of M; a typical dark matter density 
Pdm ~ lOOM /pc 3 corresponds to p DM ~ 10~ 26 A/~ 2 to 
10~ 30 A/" 2 for black hole seed masses of 10 6 and 10 4 M Q re- 
spectively. However, due to numerical restrictions in the use 
of ultralow fluid densities -as corresponds to dark matter- in 
geometrical units when solving the relativistic Euler equa- 
tions; we study the accretion of fluid densities numerically 
tractable and then extrapolate our results to the density 
corresponding to dark matter and extract our conclusions. 




Figure 1. In this figure we present the accretion mass rate and 
the accreted mass M acc for a pressure-less fluid for various values 
of the initial coordinate velocity v r = 0.0, —0.1, —0.2, —0.3. In all 
the cases, the initial density used is p = 10~ 8 . The different types 
of line indicate that the accretion is being measured at various 
spherical surfaces, and the fact that they behave similarly shows 
that the trend of the accretion rate is consistent at the various 
detectors. 



4 RESULTS 

4.1 Preparing the initial conditions 

In order to solve Euler equations (|10[) and (|11|) we have to 
provide initial data at t = 0, for p and v r that for each value 
of T would imply different values for p. In principle there 
is no prescription for p and v T which would be arbitrary. 
Since we are interested in searching steady state scenarios 
we are inclined to set rather unbiased out of equilibrium 
initial conditions in order to see whether or not the system 
relaxes and approaches a sort of late-time attractor solution. 

We have found that an initial constant density profile 
corresponds to a non-steady state condition, and instead 
produces an immediate dynamical behavior. We also choose 
different values of v r , because we want to explore the space 
parameter in this direction too. We have found that the ini- 
tial radial velocity is only restricted in the sense that for big 
initial inward velocities the fluid may reach ultra relativis- 
tic speeds when approaching the event horizon. Finally, the 
pressure is obtained from the equation of state. 

An initial parameter for which there is no prescription 
is that of the internal energy of the gas e in the equation of 
state. Its initial value determines the pressure at the initial 



time with the property that when e approaches zero the 
initial conditions corresponds to the p = case. 

4.2 Case p = 

In Fig. [T] we show the numerical results of the mass accre- 
tion rate and the accreted mass for different values of the 
initial inward velocity with an initial constant density profile 
p = 10 -8 ; the accretion rate is measured at three different 
surfaces located at r = 2M, 14M, 29M in order to make 
sure the trend of the accretion rate is not an artifact due to 
backscattering near the horizon. 

In all cases the mass accretion rate increases rapidly 
and never during our evolutions enters in a steady state. In 
order to measure how fast the accreted mass grows we fit 
the accreted mass with an anzats M = At B from our nu- 
merical results. We also extrapolate the parameters A and 
B to the cases of two reasonable density values correspond- 
ing to dark matter, and black hole seeds of 10 4 and 10 6 
solar masses. We summarize these results in Table [T] where 
we show the values of A and B that better fit our numeri- 
cal results, where the last two columns are an extrapolation 
corresponding to densities of lOOM© /pc 3 around black holes 
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Table 2. In this table, we show the mass accreted by the a black 
hole since lOGys ago to the present for two black hole seeds. The 
initial density profile of the dark matter is p = lOOM0pc — 3 , which 
in geometric units, for black hole seed masses = 10 4 Mq 

and M( U ) = 10 6 A/ Q , correspond to p = 1.17 X 10~ 30 and p = 
1.17 X 10~ 26 in geometric units respectively. The values of the 
mass accreted for these densities are taken from the extrapolation 
in Table [TJ 







seed M (i) = 1O 4 M 


seed M (ii) = lO 6 Af 
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3.04 x 10 7 M (i) 


3.82 x 10 7 M (ii) 


v r 


= -0.1 
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2.07 x 10 16 M (ii) 
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3.28 x 10 13 M (ii) 
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= -0.3 


1.96 x 10 13 M (i) 
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of 10 4 and 10 6 M Q . The domain we use in these numerical 
experiments is M < r < 501Af and made sure the running 
time t — 5000A/ is adequate for this grid size and our code 
converges. 

In order to explore the possibility that a pressure less 
fluid can be the dark matter we consider two particular ex- 
periments: 

i) we choose the system to start accreting lOGys ago, with 
an initial dark matter density of IOOA/q jpc' onto a black 
hole with mass 10 4 Mq, 

ii) and the second system starts accreting also lOGys ago 
onto a black hole seed of 1O 6 M0, 

for which the fitting parameters are available in the last two 
columns of TableQ] The results for these two experiments are 
presented in Table [2] where it is shown that the black holes 
should have accreted enough dark matter as to achieve a cur- 
rent mass of 10 11 Mq in the most conservative cases or even 
10 22 Mq in the worst case. We remind the reader that given 
that we only consider radial accretion, these results work 
fine as upper bounds for more general symmetries involving 
non-trivial impact parameters of the volume elements of the 
dark fluid. 

It can be noticed that the initial velocities explored 
range from zero up to about one third the speed of light, 
however we found similar results for even higher speeds 
nearly half the speed of light. 



4.3 Case p 

In this case, the parameter space expands because now Y 
is a new parameter and so is the internal energy e at initial 
time. In order to explore a wide parameter space we consider 
the following cases: 

(a) Again choose various values of the initial density p = 
10" 4 , 10" 6 , 10" 8 , 10" 10 . 

(b) We only study two values of the inward velocity, v r = 
-0.1, -0.2. 

(c) We choose two values of the adiabatic constant T = 
1.1, 1.2, so that we estimate self-interacting matter. 

(d) And two values for the initial internal energy e = 
0.5, 1.0, because unfortunately there is no a priori pre- 
scription for this parameter, unless for instance a polytropic 




Figure 3. In this figure we present the accretion rate of an 
ideal gas, for two values of the initial specific internal energy 
e = 0.5, 1.0 and parameters: initial density 10~ 4 , initial velocity 
v r = —0.1 and adiabatic constant T = 1.1. The role played by 
the internal energy at initial time consists in providing a value to 
the initial pressure, that is, when e tends to zero the value of the 
pressure tends to zero too. 



equation of state is assumed initially, in which case the poly- 
tropic constant would be the free parameter. 

In all the combinations of these parameters we found 
that the accretion process reaches a stationary accretion 
mass rate. We show in Fig. [2] a representative sample of 
our results indicating a linear accreted mass in time after 
a transient lapse, where the initial relax toward a sort of 
stationary late-time attractor solution. For completeness in 
Fig. [3] we show as an example the mass accretion rate with 
two different initial specific internal energies. This shows the 
role played by the initial internal energy, that is, in the limit 
of e approaching zero as can be seen in (|12[) . the pressure 
tends to zero too and one recovers the unstable case. 

Our results on the accretion mass rate and accreted 
mass are summarized in Table [3] where, as in the pressure- 
less case, we also extrapolate our fitting results to the ultra 
low density corresponding to dark matter. The domain we 
use in these numerical experiments is M < r < 1001A/ 
and made sure the running time t = 3000M is adequate for 
this grid size and our code converges. We then use these ex- 
trapolated parameters to estimate the accreted mass during 
lOGyr, which can be found in Table [4] 

From our results it call the attention that approach- 
ing a stationary accretion process is very sensitive to the 
value of F, as can be seen in the Tables [2] and [4] where 
the resulting estimates in the accreted mass can change by 
about fifteen orders of magnitude when changing T from 
1 to 1.1, which corresponds to the system of dust and our 
ideal gas dark matter model with pressure respectively. In 
order to show that this is not an artifact of our numerical 
approach, we show in Fig. [4] that the accretion mass rate 
changes monotonically with F for two configurations chosen 
from the collection of cases we experienced with. 

Long time attractor behavior of the solutions. As an 
example of how the state variables approach a stationary 
regime in the case of non zero pressure, we show in Fig. 
[5] snapshots of p, p and v r . It can be noticed that the in- 
tial data corresponding to constant density, pressure and ve- 
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Figure 2. In this figure we present the accretion rate of an ideal gas, for two values of the initial coordinate velocity v r and two values of 
the adiabatic constant T. In all the cases, the initial density and the specific internal energy used are p = 10 -4 and t = 0.5 respectively. 
The different type of line indicates that the accretion is being measured at various spherical surfaces, and it is remarkable that the 
accretion mass rate approaches the same value, as happens in the p = case, however, in the present case, due to the non-zero pressure 
the result is non-trivial because outward and inward flow are expected to happen. 
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Table 1. In this table we show the fit parameters describing the mass growth function for a wide set of initial density and inward velocity 
of the test fluid for the pressure-less case. We also show the extrapolation value corresponding to the dark matter density p = 10 -26 and 



10- 



in geometric units, that corresponds two particular cases we analyze below. 



utud mass (Dust Fluid; 



Table 3. In this table we show the parameters fitting the accreted mass for the initial internal energy e = 0.5, 1.0 and the adiabatic 
constant V = 1.1, 1.2. We also show the extrapolation value corresponding to the dark matter density p = 10~ 26 and p = 10 -30 in 
geometric units, that correspond to two particular cases we analyze below. The anzat we use to fit the mass accreted is M acc = At B + C. 
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locity, after a transient lapse approach a time-independent 
regime. It can also be observed that the pressure is bigger 
near the horizon -because it is linear with the rest mass 
density- which prevents the fluid from falling quickly into 
the hole, an effect that does not happen in the p = case 
and we consider to be the responsible of slowing down the 
accretion rates. 



4.4 Test of the code 

In order to validate our numerical results we present a very 
stringent test which is related to the convergence of our nu- 
merical solution. Since we do not have an exact solution of 
the equations, we perform a Cauchy type of self-convergence 
test on the state variables. As a sample of such test we show 
in Fig. [6] the self-convergence of the rest mass density for 
two of our runs. We find that the order of self-convergence 
is one, which is consistent with the order of approximation 
(first order) of our variable reconstruction at the cell inter- 
faces of the Riemann solver. 



5 CONCLUSIONS 

We have numerically studied the spherical accretion of a 
perfect fluid with an ideal gas equation of state on a fixed 
Schwarzschild background space-time in penetrating coor- 
dinates and support our results with self-convergence tests. 



Our procedure consists in the numerical solution of the time- 
dependent relativistic Euler equations on the black hole 
space-time. As a first approximation, due to well known 
numerical difficulties, we have extrapolated our results to 
the physical values corresponding to the accretion of dark 
matter onto a supermassive black hole. We have explored a 
wide range of parameters related to density, pressure, inter- 
nal energy, adiabatic constant and inward initial velocity of 
the fluid. 

For the parameter space explored, we have found that 
when the pressure of the dark matter is zero, within the 
window of our simulations there is no steady state accretion 
process, and also that the current masses of supermassive 
black holes should be orders of magnitude bigger than those 
observed if they are assumed to be fed by spherical accretion 
of dark matter. However, for the case of non-zero pressure, 
we have found in all the experiments we have analyzed that 
a steady state accretion is reached and behaves like a late- 
time attractor solution which is reached sooner for bigger 
values of the adiabatic index. 

We found that for initial black hole seed of already su- 
permassive size (1O 6 M0), the accreted mass during lOGyr 
is of the order of solar masses. Either this implies that dark 
matter may obey an equation of state with non-zero pres- 
sure, or in a conservative case, the accretion of dark mat- 
ter onto supermassive black holes contributes with a small 
fraction compared to possible accretion of baryonic mat- 
ter as suggested by the bolometric luminosity of quasars 
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Figure 4. We show the accretion mass rate for various values of 
r for the configurations corresponding to the initial parameters 
p = 10~ 4 , v r = —0.1, e = 0.5, 1. In the insets we show a zoom in 
of the accretion mass rate that illustrates the various behaviors 
at the initial transient stage that depends on T. 

(Pcirani & Frcitas 2008), where a bound of 10% of the ac- 
cretion is allowed to be dark matter. 

This spherically symmetric accretion model is a first 
step toward a full study of the SMBH plus dark matter 
system, which in a complete version should involve the in- 
troduction of tidal forces due to angular components of the 
fluid that will break the spherical symmetry, back reaction 
effects or full non-linear evolution of the black hole, transfer 
of angular momentum to the black hole, etc. 
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Figure 5. In this figure we show how after an initial transient 
the state variables approach a steady state, which provokes the 
accretion mass rate to be stationary. This case corresponds to the 
initial parameters p = 10 — 4 , v r = —0.1, T = 1.1 and e = 0.5. The 
set of gray, dashed and dotted lines correspond to snapshots at 
early times during the simulation, whereas the black bold line is 
actually a set of various snapshots superposed after the system 
enters a steady stage. The curves start from the excision boundary 
at r = M and we only show the region near the horizon where 
these quantities have structure. 
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Figure 6. In this figure we show the order of self-convergence of 
p, for two different physical cases, 1) (left panel) p = 0, p = 10~ 4 
and v r = -0.1 and, 2) (right panel) T = 1.2, p = 10~ 4 v r = -0.1, 
e = 0.5. We calculate the convergence using the L\ norm of the 
differences between the value of the density for three different 
resolutions Axi = 0.1, Aa;2 = Axi/2 and AX3 = AX2/2. Then, 
the self-convergece factor, Q, for these three resolutions is calcu- 
lated from the following cxprcsion as 

2 Q S c = AE1/AE2 where 
AEi = Li(pi — P2) and AE2 = £i(p2 — P3)- The dominant 
approximation is that of the reconstruction of variables at the 
interface cells, which is a piece wise constant and therefore first 
order accurate. This is consistent with the self-convergence of this 
plot. 
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